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ABSTRACT 

From an observational point of view, the role of magnetic fields in star formation remains unclear, and two main theoretical scenarios 
have been proposed so far to regulate the star-formation processes. The first model assumes that turbulence in star-forming clumps 
plays a crucial role, and especially that protostellar outflow-driven turbulence is crucial to support cluster-forming clumps; while the 
second scenario is based on the consideration of a magnetically-supported clump. Previous studies of the NGC 2264-C protocluster 
indicate that, in addition to thermal pressure, some extra support might efi'ectively act against the gravitational collapse of this cluster- 
forming clump. We previously showed that this extra support is not due to the numerous protostellar outflows, nor the enhanced 
turbulence in this protocluster Here we present the results of the first polarimetric campaign dedicated to quantifying the magnetic 
support at work in the NGC 2264-C clump. Our Zeeman observations of the CN(l-O) hyperfine lines provide an upper limit to 
the magnetic field strength i?ios i 0.6 mG in the protocluster (projected along the line of sight). While these results do not provide 
sufliciently tight constraints to fully quantify the magnetic support at work in NGC 2264-C, they suggest that, within the uncertainties, 
the core could be either magnetically super or sub-critical, with the former being more likely. 
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1. Introduction 



1.1. Background: The role of magnetic fields in star formation 

It is well-established that stars form by gravitational collapse 
of molecular dense cores, themselves embedded in parsec-scale 
clu mps contained in large molecular clouds (see, e.g. the review 
bv lLarsonll2003l) . However, the conditions under which these 
dense cores form, collapse, and fragment remain a matter of de- 
bate. In particular, understanding the physics regulating the star- 
formation processes is crucial to, e.g, reconciling the values of 
star formation efficiencies observed on galactic scales, which are 
significantly lower than expected from purely gravitational col- 
lapse JZuckerman & Evans .1974; .Krumholz & Tan ,20071) pre- 
dictions. 

Two main competing models have been proposed for driving 
and regulating the star-formation processes on the scales of star- 
forming clumps. 

The first scenario argues that molecular clouds are intermittent 
struc tures in an inters tellar medium dominated by turbulence 
(e.g. lElmegreenl |200Cl), and that turbulent motions prevent the 
clouds from c o llapsin g in a freefall time (Mac Low & Kl esseni 
l2004t iMversI l2000l) . In this m odel, star-forming clumps 
undergo turbulent frag mentation ( iPadoan & Nordlundl 120021 : 
IClark & Bonnelll l2005h : self-gravitating pre-stellar condensa- 
tions form as turbulence-generated density fluctuations, turbu- 
lence then dissipates rapidly, and the cores eventually collapse 



* Fits file for the CN(l-O) map is only available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http: //cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 

** Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



with little interaction with their surroundings. A somewhat dif- 
ferent p icture has been proposed in the case of clustered star for- 
mation, amura '&Lil(l2007l) argued that, owing to its short 
decay time dMac Low etaT 119981) . the interstellar turbulence 
that is initially present in a cluster-forming cloud is quickly re- 
placed by turbu lent motions generated by the numerous proto- 
stellar outflows. iLi & Nakamural (l2006l) proposed that the proto- 
stellar outflow-driven turbulence dominates for most of a proto- 
cluster's lifetime and acts to maintain the cluster-forming region 
close to overall virial equilibrium for several dynamical times, 
avoiding global free-fall collapse and reducing the local star- 
formation efficiency. 

The second model proposes that the magnetic field provides 
an efficient support against gravity and delays the star forma- 
tion (Shu et al. 1987^, by forming self-gravitating dense clouds 
that are magnetically supported (e.g. , Mouschovias & Spitze^ 
119761: ICiolek & Basull200lt iHennebelle & Tevssierii2008,) '. The 
crucial parameter characterizing magnetic support of a dense 
cloud against its self-gravit y is the ratio of the mass to mag- 
netic flux M/<t>B. Mouschovias & Spitze^ ( 1 19761) calculated the 
critical value of this mass-to-flux ratio using the virial theorem 
coupled to numerical calculations of a bidimensional cloud, and 
found that (M/Ob)^/ ~ 0.13 x G"'^^, where G is the gravita- 
tional constant. A cloud that has a mass-to-flux ratio larger than 
this critical value is said to be magnetically supercritical: it can- 
not be supported by magnetic pressure, and collapse can pro- 
ceed on clump-sized scales, gravitationally bound dense cores 
eventually being formed. The cores progressively evolve towards 
higher degrees of central concentrations as the magnetic support 
is progressively lost through ambipolar diffusion. Otherwise, the 
cloud is subcritical: gravitational collapse is impeded because 
magnetic fields support the cloud against self-gravity. It is how- 
ever possible that supersonic turbulence accelerates star forma- 
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tion in clouds that are initial ly subcritical, through enhanced am- 
bipolar diffusion in shocks (iNakamura & Lill2005h . 

One of the reasons for the ongoing debate about the re- 
spective roles of turbulence and magnetic fields in regulat- 
ing star formation is the difficulty in accurately measuring 
the magnetic field in the dense environments typical of star- 
forming clumps. On the observational side, only a handful of 
magnetic field strength estimates have been reported at den- 
sities i 10'' cm"^. While low-density HI structures (such as 
Giant Molecular Clouds) are found to be globally magneti- 
cally subcritical CHe iles & Crutcher 2005 ). most observations 
towards dense molecular cloud cores on sub-pc scales ob- 
tained so far suggest t hat th ey are approximately critical or 
supercritical ( ICrutched l2009l) . Typically, magnetic fields de- 
tected in t hese dense molecular clouds have strengths 0.1 
1.5 mCi (ICrutcher et al l 1 19991 l2004t iFalgarone et all l2008t 



iTroland & Crutchej|2008|- iHever & Bruntll2012h . We note, how- 
ever, that some recent work based on the statistical analysis of 
the linear polarization suggests that some cloud cores (or re- 
gions of cloud co res) may be subc r itical, such as the b owl of 
the Pipe Nebula dAlves etalJ 120081: iFranco et al.]l2010l) or the 
massive star-forr ning reg ion W5 1 , on scales ~l-2pc ( Tang et aP 
l2009tlKochetani20121 



1.2. Previous studies of fhe NGC 2264-C protocluster 

The NGC 2264-C cluster-forming clump is located in the 
Mon OBI giant m olecular cloud complex at d~ 800-900 pc 
( Baxter et al.l 120091) . Performing 30m and PdBI observations, 
iPeretto et al.l ( |2006|) carried out the first comprehensive mil- 



limeter continuum/line study of the protocluster. Their 1.2 mm 
continuum mosaic of NGC 2264-C resolved the internal struc- 
ture of the region, uncovering a total of 13 compact prestel- 
lar/protostellar cores. Moreover, their line observations, com- 
bined with radiative transfer modeling, established the presence 
of large-scale collapse motions, converging onto the most mas- 
sive core (C-MM3, near the center of NGC 2264-C ). Detailed 
comparison of these observations with numerical smoothed- 
particle hydrodynamics (SPH) simulations of the evolution of a 
IOOO-M0 Jeans-unstable isothermal clump dPeretto et al.ll2007h 
suggested that NGC2264-C is an elongated clump collaps- 
ing/fragmenting along its long axis, and at a very early stage 
of global clump collapse ( i lO*" yr after the start of dynami- 
cal contraction). A significant shortcoming of their SPH simu- 
lations, however, is that they only produce the observed level of 
clump fragmentation when the total mass of dense (> 10"* cm"^) 
gas in the model is a factor of ~ 10 lower than in the actual 
NGC 2264-C clump. This suggests the existence of forces that 
are not included in the hydrodynamical simulations and act as 
an extra support against gravity in NGC 2264-C, such as some 
support provided by feedback from protostellar outflows or mag- 
netic fields. 

To test the hypothesis o f outflow-generated s uppor t sug- 
gested by iLi & Nakamural (l2006h . iMaurv et al.l (l2009h used 
the heterodyne array HERA on the 30-m telescope to con- 
duct a search for protostellar outflows, by mapping '^CO(2-l), 
'^CO(2-l), and C"^0(2-l) over the protocluster. They found a 
total of eleven outflow lobes originating from the dense proto- 
stellar cores forming in the protocluster They carried out a quan- 
titative study of the momentum flux injected by these outflows 
into the protocluster, and concluded that the network of outflows 
at work in NGC 2264-C fails to efficiently support the protoclus- 
ter against global collapse. 



The NGC 2264-C protocluster is therefore a well-studied 
star-forming clump, showing relatively simple global-collapse 
motions. To improve our understanding of this object but also 
use it as a case study for testing theories of clustered star for- 
mation, the magnetic field strength needs to be quantified in this 
dense clump. This should allow us to determine whether some 
magnetic support is responsible for the observed properties of 
the cluster-forming clump, and test the aforementioned theory 
of magnetized clouds. In this paper, we present the results of the 
first polarimetric campaign dedicated to estimating the magnetic 
field strength in the NGC 2264-C cluster-forming clump. 



1.3. CN Zeeman effect 

If a medium is permeated by a magnetic field B, the energy lev- 
els of the molecules and atoms split, owing to the interaction be- 
tween the magnetic field and the magnetic dipole moment asso- 
ciated with their orbital angular momentum. The splitting of the 
energy levels then results in the splitting of the spectral lines into 
several separate polarized co mponents : this e ffect, which was 
observed for the first time by IZeemanI ( Il897h . is the so-called 
Zeeman effect. 

The Zeeman effect remains the only direct method for 
measuring magnetic fi eld strengths in dense molecular clouds 
(ICrutcherl l2007l l2009l) . but the range of molecules that are at 
the same time abundant, paramagnetic, and accessible at mil- 
limeter wavelengths is rather small. As a tracer of high-density 
gas, the CN thermal lines probe the densest regions in molecu- 
lar clouds. Moreover, the N = 1^ transition of CN consists 
of nine hyperfine components (hereafter HFS), of which seven 
are strong components sensitive to the Zeeman effect. These 
HF S lines have different Zee man splitting factors (see Table[T] 
and [Turner & Gammonll 19751 for further references). This prop- 
erty of the CN(l-O) emission lines is crucial for performing ac- 
curate Zeeman measurements probing the magnetic field, since 
the magnetically less-sensitive HFS (components Z2, Z3, and 
Z6 in Table [U can be used to c orrect for sidelobes a nd in- 
strumental polarization effects (see ICrutcher et al1ll996l for the 
detailed procedure). Since these instrumental effects can some- 
times be stronger than the Zee man signal produced by the mag- 
netic field in the source itself dForbrich et al.ll2008l) . the ability 
to make these corrections is crucial to distinguishing false detec- 
tions from Zeeman signals that are truly magnetically-induced 



Table 1. CN(l-O) strong hyperfine lines 



Line 


Frequency 


Relative 


Zeeman splitting 


|RI X Z| " 




(GHz) 


intensity RI 


factor Z (Hz/^iG) 




Zl 


113.144 


8 


2.18 


17.4 


Z2 


113.171 


8 


-0.31 


2.5 


Z3 


113.191 


10 


0.62 


6.2 


Z4 


113.488 


10 


2.18 


21.8 


Z5 


113.491 


27 


0.56 


15.1 


Z6 


113.500 


8 


0.62 


5.0 


Z7 


113.509 


8 


1.62 


13.0 



" |RIxZ| is the relative sensitivity to the magnetic field along the line of 

sight Sfoj. 
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RA (J2000) 

Fig. 1. Integrated intensity map of the CN (A^ = 1 - 0, 7 = 3/2 - 1/2) hyperfine components, obtained with EMIR in May 2009. 
The black contour represents the 3-cr integrated intensity level (7.3 K .km.s ') in the map . The crosses show the location of the 
millimeter dense cores mapped with MAMBO on the IRAM-30m by iPeretto et al.l (l2006l) . The white d ashed contour show s the 
value of column density ~ IxlO^"* cm inferred from the MAMBO and P-ArTeMiS continuum maps dMaurv et al.ll2009l) . The 
circle shows the location "Z" where we carried out further polarimetric observations to probe the Zeeman effect. 



2. Observations and data analysis 

2.1. CN mapping of the NGC 2264-C protocluster 

iPeretto et al.l (l2006h showed that the typical densities of the 
dense material in NGC 2264-C range from a few 10^ cm""* to 
a few 10^ cm"^ With its critical density n{H2) ~ 10^ cm ^ CN 
is therefore an appropriate candidate molecule for conducting 
Zeeman observations in this region. 

In preparation for future CN Zeeman observations, we first 
carried out an extensive mapping of the cluster-forming clump 
in the CN(l-O) line, to find positions where the line strength and 
profile would allow some Zeeman measurements. These prelimi- 
nary observations were carried out in May 2009 with the IRAM- 
30m telescope (FWHM of 23" at 113 GHz). We mapped a region 
of 5' X 3' centered on C-MM3 (see Fig.[T]), using the on-the-fly 
mode. The Eight Mixer Receiver (EMIR) frontends were tuned 
to 113.3 GHz, thus allowing us to observe simultaneously both 
the first (Zl, Z2, Z3) and the second (Z4 to Z7) group of CN( 1- 
0) hyperfine components. The use of the VErsatile SPectrometer 
Assembly (VESPA) backends resulted in high spectral resolu- 
tion (80 kHz i.e 0.1 km/s) maps for each group of HFS compo- 
nents. We used simultaneously a third VESPA backend to cover 
a wide spectral window of 640 MHz, with a lower spectral res- 
olution of ~0.8km.s"'. With 12 hours of total observing time, 
we obtained root mean square (rms) noises in the final maps 
of -200 mK per 0.8 km.s^'-channel. Fig. [1] shows the map of 
the CN(l-O) emission, integrated over the second group of HFS 
components (between 113.48 GHz and 113.52 GHz, transitions 
with J - 3/2 - 1/2). The CN emission is strong towards the 
protocluster, and closely follows the dust continuum emission. 

The spectra obtained towards the high density central region 
of the protocluster show double-peaked profiles, which could 



be due either to self-absorption, infall motions or the presence 
of two velocity components. To be able to distinguish between 
these three scenarios, we would need to perform observations in 
another optically thin, single-peak ed transition of CN. We n ote, 
however, that earlier observations ( IPeretto et al.ll2006l l2007h us- 
ing other molecular tracers of dense gas have discovered the 
presence of large-scale infall motions in the protocluster. This 
strongly suggests that the double peak spectra are due to infall 
motions, but does not exclude other scenarios. 

In the external parts of the protocluster, the CN emission is 
well-represented by single-peaked Gaussian line profiles, which 
are indicative of an optically thin emission. Moreover, the spec- 
tra obtained in these external parts (see Fig. [2]) show that al- 
though the measured relative intensities of the seven HFS com- 
ponents vary slightly from their optically thin LTE ratios (re- 
ported in Table[Tli, the variations are not large, thus suggesting 
that the CN emission is optically thin towards the outer parts of 
the clump. Finally, we stress that the components "0" and "8" 
shown in Fig. |2] are respectively the CN hyperfine components 
with quantum numbers (J=l/2-l/2, F= 1/2-1/2) and (J=3/2-l/2, 
F= 1/2-3/2), not used in the Zeeman analysis because their inten- 
sities are too weak. 

2.2. Polarimetric observations 

Following iFalgarone et al.l (l2008l) . we estimate that CN(1- 
0) lines stronger than ~1K enable CN Zeeman observations 
sensitive to fields Bios~0.2 mG within a moderate amount 
( i 40 hours) of telescope time, u sing the XPOL polarimeter on 
the IRAM-30m telescope (iThum e t al. 2008). The CN emission 
has to be close to optically thin, however, to avoid confusing the 
effects of optical thickness in the beam with a Zeeman signal. 
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Fig. 2. CN(l-O) spectra from the "Z" position used for Zeeman observations, highlighted in the CN(l-O) integrated intensity map 
here-above, (a) Broad-band spectrum obtained at low spectral resolution (0.8 km. s"' channel), covering the complete extent of 
all CN hyperfine components, (b) VESPA spectrum obtained towards the CN(l-O) triplet (J=l/2-l/2) at high spectral resolution 
(O.lkm.s ' per channel originally, here smoothed to a 0.2km.s ' channel), towards the Z position. A fit of the hyperfine structure 
(red dashed spectrum) indicates that the systemic velocity is Vlsr - 8.7 + 0.4 km/s. (c) A VESPA spectrum obtained towards 
the CN(l-O) quadruplet (J=3/2-l/2) at high spectral resolution (O.lkm.s"' per channel originally, here smoothed to a 0.2km.s"' 
channel), towards the Z position. 



The "Z" position (ajzooo = 06Hl'"10.0\ dnooo = 09°28'50.00") 
highlighted in Fig.[T]is located near the center of the protocluster, 
and shows a strong non self-absorbed CN(l-O) line {T^ ~ 2 K 
in the strongest hyperfine component, see Fig.|2};). It was there- 
fore chosen to perform deep polarimetric observations probing 
the Zeeman effect in the CN lines. 

The polarimetric observations with XPOL consist in using 
two orthogonally polarized 3-mm heterodyne receivers to ob- 
tain both the auto-correlation of each receiver (horizontal H and 
vertical V spectra) and the polarization cross-products (real and 
imaginary parts of the cross-correlation spectra). We used the 
3mm receivers of EMIR associated with the VESPA backends in 
the polarimetry mode to carry out position-switch observations 
of the Z position, using an off position located at (-160", -30"), 
checked to be clean of CN(l-O) emission. The VESPA correlator 
was split into two separate parts, both with 80 kHz channel spac- 
ings and tuned to obtain simultaneous spectra of the two groups 
of CN(l-O) HFS components. The spectra were first converted to 
a temperature scale by applying the standard calibration proce- 
dure for spectral observations (ambient load, cold load, and sky 
calibration). The correction for the phase difference between the 
receivers was done by using a polarization g rid to observe the 
cold load during the calibration procedure (see lThum et al 1120081 
for further details), resulting in a phase accuracy of ^ 0. 1° rms. 
We checked the phase drift during our observations, estimating 
it to be less than 0.5° per hour Finally, we fitted and subtracted 
the baselines with a polynomial of order 1, simultaneously for 
the two VESPA parts tuned to the two groups of hyperfine com- 
ponents. 

To obtain significant constraints on the magnetic support 
at work in NGC 2264-C, we aimed to achieve a sensitivity 
o'rmsiTD - 0.5 mK in the Stokes V spectrum (imaginary part 
of the cross-correlation between the two polarizations, see §4), 
and were granted a total of 80 hours of IRAM/30-m telescope 



time to perform polarimetric observations at the Z position. The 
observations took place over one week at the end of April 2010, 
in mixed weather conditions (water vapor ranging from 4 mm to 
12 mm and system temperatures from 200 K to 700 K) that did 
not allow us to reach the desired rms, but to obtain nevertheless 
good quality data suitable for a first meaningful analysis, that we 
present here. 

3. Results: upper-limit on Bi^, in NGC 2264-C 

Our observations provide us with the Stokes parameters. The 
Stokes I spectrum is the sum of the horizontal and verti- 
cal auto-correlation spectra / = H + V, while the circular- 
polarization component Stokes V is the imaginary part of the 
cross-correlation spectrum. In the case of a magnetic field Bjos 
that is stronger than our detection limit, a Zeeman effect should 
be detected as an S-shaped signal of the form V{v) - Bios x 
Z{dl/dv) in the Stokes V spectra of the hyperfine CN lines that 
are sensitive to magnetic fields (components Zl, Z4, Z5, and 
Z7). The insensitive components Z2, Z3, and Z6 should feature 
only noise if no significant instrumental effects are present. 

The cumulative Stokes I and Stokes V spectra obtained for 
the Zeeman sensitive CN(l-O) hyperfine components are shown 
in Fig. [3^, b, and c. To improve the signal-to-noise ratio and 
therefore the detectability of the magnetic field, we computed 
the average (weighted by the sensitivity to the Zeeman effect) 
Stokes I and Stokes V spectra derived from components Zl, Z4, 
Z5, and Z7. The average spectra are shown in Fig.[3}l- Our obser- 
vations do not detect any Zeeman signature: the average Stokes 
V spectrum contains only noise, and no significant S-shaped sig- 
nal is detected. We note that the average Stokes V spectrum com- 
piled from the insensitive components Z2, Z3, and Z6 can be 
well-fitted by a single Gaussian noise component (see Fig. [3^), 
showing therefore that the instrumental effects are neglictible at 
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Fig. 3. (a) to (c) panels show the Stokes I (top panel) and Stokes V (bottom panel) cumulative spectra obtained for each of the 
Zeeman-sensitive CN(l-O) hyperfine components: component Zl in (a), components Z4 and Z5 in (b), and component Z7 in (c). 
Overplotted on each Stokes V spectrum is the measurement of dl/dv computed from the Stokes I spectrum, and scaled to the upper- 
limit magnetic field Bio, mG that can fit to each individual hyperfine Stokes V spectrum. The typical rms noise level in these 
spectra is ~2.5 mK. The panel (d) shows the average Stokes I (top) and Stokes V (bottom) spectra, using the four components that 
are most sensitive to the Zeeman effect (Zl, Z4, Z5, and Z7) weighted by their respective Zeeman sensitivity. The dashed colored 
lines overplotted in the lower panel show the dl/dv computed from the average Stokes I spectrum and scaled to match the 3-cr rms 
noise level in the Stokes V spectrum. The rms noise level in the average Stokes V spectrum is o-„t^s{T'^)~13 mK, leading to an upper 
limit on the magnetic field strength along the line of sight of Z?io.s i 0.6 mG. The rightmost panel (e) shows the average Stokes I 
(top) and Stokes V (bottom) spectra, using the three components least sensitive to the Zeeman effect (Z2, Z3 and Z6) weighted 
by their respective Zeeman sensitivity. These components are mainly sensitive to the effects of instrumental polarization, and their 
spectra do not contain any significant feature. 



the sensitivity reached by our observations. To illustrate the sig- 
nal that would be observed in the Stokes V spectrum if a suffi- 
ciently strong magnetic field were present, we have over-plotted, 
in Fig. 3d, the dl/dv signal computed from the average Stokes 
I spectrum and scaled for magnetic field strengths at our sen- 
sitivity limit B|o,s ~0.6-0.8mG (respectively, levels of 3-cr and 
4-0-). 

4. Discussion 

It is usual to define the ratio of the observed to critical mass-to- 
flux ratios /i = [M/<l)]ob,s/[A^/'5]crit as the characteristic quantity 
measuring the stability of magnetized clouds: values of // above 
1 correspond to cases where the magnetic field cannot sup- 
port the cloud against its self-gravity and gravitational collapse 
will p roceed. The theoretical study of iNakano & Nakamural 
(Il978h determined ju =7.6 xlO""'* N{H2)I\B\, where N{H2) is 
in cm"^ and \B\ is the magnetic field strength expressed in 
mG. Considering the mean column density of NGC 2264- 
C (A^Ho~l-8xlO^-'cm"-), the equipartition of the magnetic and 



gravitational energies is reached for a total magnetic field 
strength |B| ~1.4mG: such a magnetic field (or stronger) would 
be able to prevent the global gravitational collapse of the whole 
NGC 2264-C clump. However, a magnetic field strength of half 
this critical value (// - 2) is enough to provide significant sup- 
port, slow down the gravitational collapse , and might explain the 
motions observed bv lPeretto et alj (l2007h . 

To efficiently constrain the magnetic support in NGC 2264- 
C, we initially aimed to probe mass-to-flux ratios ju i 2, 
corresponding to total magnetic field strengths B i 0.6 mG. 
Assuming a statistical mean inclination of the magnetic vector 
~57° with the line of sight (or ap plying the mean statistica l ge- 
ometrical correction proposed bv iHeiles & Crutcheill2005l) . the 
line-of-sight component probed by the Zeeman effect B\o^ is half 
the value of the total magnetic field vector norm \B\. Therefore, 
we aimed to reach a typical sensitivity limit 3crg,^^ i 0.3 mG, 
corresponding to sensitivities "io-^rnAT*^ ^1.5 mK in the aver- 
age Zeeman-sensitive (components Zl, Z4, Z5, and Z7) Stokes 
V spectrum. 
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Owing to the variable weather conditions during our 30-m 
observations, we could only reach 3cr„j^s{TD - 3.6 mK in the 
average Stokes V spectrum (see Fig.O, i.e. our observations pro- 
vide an upper-limit to the magnetic field strength projected along 
the line of sight of 3xcrj„s(Z?ios)~0.6 mG. C orrecting the mag- 
netic vector for statistical geometrical effects (iHeiles & Crutcherl 
I2005h . our result turns into a tentative upper-limit to the field 
strength of Z? ;^ 2x3crrn,s(Bios)~ 1-2 mG. This is marginally 
above the dividing line between the magnetically subcritical and 
supercritical cases (yu ~ 1.1), suggesting that the NGC 2264-C 
cloud might be magnetically supercritical. 

We note that, at the position chosen for our Zeeman measure- 
ments, the CN line intensity is subject to a slight gradient across 
the 23" half-power beam width (HPBW) (see Fig.[T]i, which we 
estimate to be ~l-2 K, increasing along a south-north axis. This 
implies that the instrumental Stokes V will get a sign, depending 
on the sidelobe (the positive or negative one) where the emission 
is the strongest. However, since this intensity gradient is similar 
on all velocity channels, it cannot mimic a Zeeman feature or ar- 
tificially erase one: it is therefore unlikely that the CN intensity 
gradient could affect our Zeeman measurements. In addition, the 
Z position is located at the border of the blue outflow lobe F2 and 
in the red outfl o w lob e Fl 1 traced in the '^CO emission maps 
of iMaurv et al.l (|2009|) . We note that the CN profiles obtained 
in Stokes I might partly trace some outflowing gas. This would 
explain the rather structured velocity profiles of the CN lines 
seen for example in Figs. 2 and 3 which both have a prominent 
blue wing and a fainter red shoulder seen in the strongest com- 
ponents. These velocity features do not jeopardize our Zeeman 
analysis, however, since we search for Zeeman signatures near 
the emission maximum, i.e around the systemic velocity of the 
cloud where the emission is dominated by quiescent material and 
therefore outflowing material should not significantly contami- 
nate neither the observed Stokes I nor Stokes V. 

However, we stress that several factors can lead to the non- 
detection of a Zeeman signal even in the case of strong mag- 
netic field. First, unfavorable inclination eff'ects can result in a 
very faint line-of-sight component of the magnetic field, while 
the absolute strength of the magnetic vector could be significant 
and even suff icient to affect t he clu mp dynamics on large scales. 
For example. lOstriker et al.l (1200 ll) suggested that mean line-of- 
sight magnetic field strengths may vary widely across a projected 
cloud. It is also possible that complex motions in the protoclus- 
ter act to cancel or weaken any Zeeman signal in our obser- 
vations of the CN(l-O) emission lines, obtained in a relatively 
large beam (23" HPBW, corresponding to a pr ojected d i stance 
of 0. 1 pc in the plane of the sky). For example, iLi et al. l (l20Tlh 
recently showed that strong coupling between turbulent gas and 
the magnetic field can lead to complex magnetic field morpholo- 
gies that are strongly dependent on the turbulence properties and 
the local density. 

Deeper Zeeman observations allowing us to reach weaker 
levels of polarization (down to ~Bios i 0.1 mG, similar to 
the magnetic fie ld strength obser ved in the nearby clump 
NGC 2264-D bv iKwon et al.l 1201 ih are needed to confirm our 
first results presented here. Additional observations probing the 
component of the magnetic vector in the plane of the sky, 
such as dust polarization observations in either the infrared 
(ISugitan i et al.ll20l"ol: Chapman et al.ll201 ll) or the subm illimeter 
regime jMatfliews & WilsonI 120021: iKoch et al.ll2012ah are cru- 
cially needed to definitively rule out the possibility of strong 
magnetic support in NGC 2264-C. 
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